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ABSTRACT 


Two dimensional distortions of a panel structure due to 
welding are analyzed by the finite element displacement method 
with an assumption of elastic deformation during welding 
process. Construction of stiffness matrices relevant to the 
welding problem is demonstrated for the one and the two 
dimensional cases. Computer programs are presented. 

Computations have been carried out using the one 
dimensional experimental value of unconstrained angular change 
along the welded edge and its equivalent constrained welding 
moment aS an input to the computer program. 

These computed results are compared to the particular 
experimental values. 3!°¢ results indicate that an analytical 
approach using fir <u. -lement method can predict the distor- 
tion phenomena of a complicated two dimensional structure with 


a reasonable degree of accuracy. 





Several recommendations are made concerning further 
investigations aimed at evaluating relevant value of angular 


change and equivalent constrained welding moment. 
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NOMENCLATURE 
Length of a finite element in x- and y- 
direction, respectively 


Equivalent welding moment for one 
dimensional deflection 


Equivalent welding moment along y-direction 
as in Figure 6 


Equivalent welding moment along x-direction 
as in Figure 6 


Coefficient matrix 

Matrix defined in equation (49) 
Matrix defined as Dey es 

Young's modulus of elasticity 
Matrix defined in Appendix A 


Transpose of the matrix ee F 


FF, G, H and q, 
respectively. Y 


BD, Va fia. 


Stiffness matrix for one rectangular finite 
element 


Additional stiffness matrix due to welding 
defined in equation (58b) and (58c), 
respectively 


Length of one panel structure in x- and 
y-directions, respectively 


Number of finite element in x- and y- 
directions, respectively 


Matrix defined in equation (64), (65), (66) 


and (67), respectively 


Load matrix due to welding defined in 
equations (58d) and (58e), respectively 


Thickness of plate 


Matrix defined in Appendix A 


x 
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inverse Matra xo 

Total strain energy 

Elastic strain energy 

Welding energy 

Deflection in plate 

Weight of electrode consumed, gram/cm. 
Weight for numerical integration 

Ordinally coordinates for a finite element 


Numerical integration points in normalized 
coordinates 


Constant defined as ee 
Poisson's ratio 

Free joint angular change 
Constrained angular change 


Variational notation 








dal: 
I INTRODUCTION 


A. Background 


Welding is used extensively in the fabrication of many 
structures, including ships, eee ce buildings, pressure 
vessels, etc., providing many advantages over the techniques 
such as riveting, casting, and forging. However, distortion 
problems are always encountered during welding. Distortion 
often prevents yeeeeuemewent-ecr desicn tolerances, reduces 
joint strength by mismatching, and imparts initial deflection 
into structural members. 

A typical structural member is a panel structure, which 
is composed of a flat plate and longitudinal and transverse 
stiffeners fillet welded to the plate, as shown in Figure l. 
Angular changes produced along the fillet welds cause out-of- 
plane distortion of the panel. The excessive out-of-plane 
distortion reduces the buckling strength of the panel that 


(3) 


is subjected to compressive loading. Corrugation failures 
of bottom plating in some welded cargo ships are believed to 
be due primarily to the reduction of buckling strength of the 
plating with excessive initial distortion, 22716717) 

Welding distortion is very complicated phenomena 
associated with many fac“ors:. 

1. Heat conduc'siv! Ly of weldment which controls the 

ais er neon ee temperature gradient of welded metal. 


2. Thermal expansion coefficient which controls the 


expansion and shrinkage of welded metal. 


Ss 
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Ihe: 

3. Moduli of elasticity which governs ethe siete, 
against deformation. 

4, Yielding strength which governs the size of the 
plastic zone near the welds. 

5. Degree of constraint which governs the freedom of 
movement during weld deformation. 

6. Amount of, heat input used which governs the size of 
the melted zone. 


io AMOuUNnt of weldment used. 


Therefore, it is extremely difficult to analyze the 
welding deformation phenomena by a pure analytical way. Yet, 
only one dimensional deformation of panel structures has been 
analyzed with the assumptions of elastic deformation during 


(1) 


Peiaing process. Obviously, welding deformations are 
neither pure elastic nor pure plastic phenomena, but the 
plastic zone near the weldments is so small compared to panel 
structural size that it is believed to be closer to the 
elastic behavior. 

Therefore, the analytical approach of the one dimensional 
deformation prediction can be treated as an elsstic deforma- 
tion by considering moduli of elasticity, degree of constraint, 


and amount of weldment used. 
B. Previous Investigations 


Masubuchi, et. al. 113) have investigated the one 
dimensional plate deformation during welding with a minimum 


strain energy concept using elastic bending theory. 
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Figure 2 Distorticon-Caused by Angular Change in Two Types 
of Fille’. Welded Structures 
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If a fillet joint is free from external constraint, the 
joint simply bends to a polynomical form having a knuckle 
at the weld as shown in Figure 2a. However, if the joint is 
constrained by some means, a different type of distortion is 
produced. For example, if the stiffeners are welded to a 
rigid beam, as shown in Figure 2b, the angular changes at 
fillet welds cause wavy or arc-form distortion. 

In the simplest case in which sizes of all welds are the 
same, the distortions of all spans are equal and the distor- 
tion, W, which can be expressed as follows: 

y= 174 = (%/2 = 1/2)7) - 6 (1) 
where 6 = deflection of plate 


Q angular change at a fillet weld with constraint 


il 


~ 
i 


length of span 


The amount of angular change, 9, in constrained 
structures is smaller than that in a free joint, 86° This 
indicates that a certain amount of energy is necessary to 
decrease the angular change from a to 8. If the necessary 
energy 1S represented by Uae it may be given by the following 


equation. 


8578 aU 
a STS oy d (8-6) (2) 


On the other hand, the strain energy stored in the 
constrained plate, Da can be expressed uSing the elastic 
beam theory. 

i ee 








16 
where E' = E/l - oe 


fe Eo 


Since U,, increases but on the contrary ae decreases as 
the constrained angle 6 increases, the condition for equili- 
brium of this system requires that the total strain energy 
U, = Ce + Le Should be minimum. Furthermore, for the sim- 
plicity of the problem, the ratio of incremental welding 
energy change to angular change is assumed to be linear as in 


the following equation. 


qu, 


Ao 6) = C(@ 8) (4) 
238: 
where C = a coefficient determined by the weight of the 
deposited metal and by the welding procedures. 


From equations (2) and (4), welding energy per unit 


length of width can be expressed in terms of C, 86, and OO: 


Bric ote 
= > (8 | 0) (5) 


Accordingly, the condition of equilibrium is as follows: 


dU 


Ww _ Le 
ag = C(O am = (6) 


From equation (6), relation between 8 and 86 can be 


expressed, 
6 


O 
eae (7) 
ie aE /c 





The values of C can be determined by using relationship 
obtained by experimentally measuring on and "6s itip died 


formula for the calculation of C has been proposed by 


MS 
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Figure 3 Variation of Angular Change of Free Fillet Welds, 
85, As a Function of Plate Thickness, t, and 
Weight of Electrode Consumed per Weld Length, W. 








Alpe: 


Peeabuchi. ‘>? 


@ Bet ems kg-mm/mm radian (8) 


where W = weight of electrode consumed, gram/cm. 


I 


2 thickness of plate, mm. 


The experimental results of determining the unconstrained 
angular change as a function of plate thickness, t(mm), and 
weight of electrode consumed per weld length, W(gr/cm), have 


been given by Hirai, et. al. \? 


aS in Figure 3. 
Therefore, for a given plate geometry and welding con 
condition, the deflection of the one dimensional case can 


easily be estimated by using equations (1), (7), (8), 


ena Figure 3. 
C. Aim and Purpose of Present Studies 


As stated previously, welding deformation has been 
analyzed only for the one dimensional case, but in reality 
most of the structural members are encountered in two 
dimensional deformation, in which the deformation, w, varies 
along the x- and y-directions as well as the unconstrained 
angular change, 6, and equivalent constrained welding moment, 
C. With the idea of previous investigations which state that 
the total strain energy in the plate during welding deforma- 
tion becomes minimum, it is afaeiene: WS extend the one 
dimensional deformation analysis to the two dimensional case. 
However, the simple mathematical solution using plate bending 


Equation is rather difficult due to the uncertainty of the 





eS 
equivalent load and the boundary conditions relevant to the 
welding deformation problem. 

The difficulties can be easily overcome by formulating 
the problem into a finite element model because it can be 
self-adjusted to take care Of the equivalent ‘Were and eguiva- 
We sittfness matrix due to welding. Furthermore, in this 
way the relation of deformation phenomena to the constrained 
angular change, 9, and the equivalent constrained welding 
moment C can be easily seen. 

Therefore, it is intended to handle the two dimensional 


deformation with the finite element displacement technique. 


\ 
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If FORMULATION OF FINITE ELEMENT EQUATION 


A. One Dimensional Case 

In this study, computer programs are programmed only 
for the two dimensional finite element analysis, but the 
one dimensional finite element formulations are presented to 
visualize the inside of the complex matrix assembling and 
computational procedure; furthermore, to compare the results 
to that of the previous investigators. 

For the one dimensional case, as shown in Figure 2b, if 
deflections are elastic and small, the deflection can be cal- 
culated using elastic theory when ae welding equivalent loads 
and the boundary conditions are prescribed. Even though the 
deformations are assumed to be elastic for the welded struc- 
tural case, the loading and boundary conditions to be applied 
for this problem are not clear. | 


As introduced by Masubuchi, et. al. ‘2) 


, if the concept 

of equivalent constrained welding moment, C(lb-in/in rad.) 

is being used, the uncertainty of loading conditions may be 
overcome; for the boundary condition it may be safe to assume 
to be simply supported for which only the angular change along 
the welded edge can be allowed. Furthermore, EO Simp leEGrtiry, 
of the problem, the deflections are considered only for the 
case of one span length of a simply supported beam to elimi- 


meee the complexity of interaction due to the sStatically 


indeterminate structural effects, as in Figure 4a. 








ail 





‘a. Simplified one-span beam element 





b. Local coordinate and nodal points for a 
finite beam element 


Figure 4 One-span Beam Element and the Finite 
Element Local Coordinate 
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For the one dimensional plane stress, the strain energy 
can be expressed simply by: 


Os = 1/2 ffs €,,9,,dxdydz (9) 


Using the stress-strain and strain-curvature relation for the 


beam element, 


= EEC (10a) 
x Xx 


2 
c= -z2 28 | (10b) 
dx | 


Therefore, the strain energy can be expressed in terms of 


deflection and rigidity of the plate: 








a Z 2 
5. = 2? | (11) 
P 2 0 Bee 
where 
q ees S 
me > OS f z*dzay = 2 
0 0 


On the other hand, the welding energy can be expressed 
in terms of ae 86, and C, as in eguation (5). For this 
problem it 1S more convinient to express the welding energy 
in terms of deflection rather than the angular changes. From 
equation (5 


Cc 2 
5 (90-8) dydx 


hs 
a 
oe = ff 
ee 


on 


the constrained angular change can be expressed by the 


equation 
8 = dw/dx (12) 


Therefore, the welding energy expressed, 





D7) 


a 2 
BEDE _ ow 
a) : (9, - 33) dx (13) 


Where the values of C and 8 A are constant which will be 
determined by the particular welding condition and the 
geometry. 

The equilibrium condition for the system requires that 
the total strain energy has to be minimized, in other words, 
the variation of the total strain energy has to be zero, 
which becomes 


SUp = 5U,, + SU, = 0 (14) 


using equations (11) and (13), equation (14) becomes 


a 1 D a at 
moe 6/22 pf (2 ax! 4+ 8 Bo f (8o- Gy)” ax|= 0 
Ab a) ' dx7 ax 


Bue Ehe variation of the constant term is zero, therefore, 


a 


a 2 dew a a 
da“ ow be a dow 
ers () (Sax + BS 5 (SY) (Bax - be s (Sax (15) 
0 0 0 
For a beam element as shown in Figure 4b, the deflection 
can be assumed to be a general cubic function of x, 
2 3 


w= cy +4 Cox + C3x + C4x (16) 


and when the nodal displacements and the nodal angles are 

the unknown variables which are to be calcualted, the coeffi- 
Cients (cy, Cor Cay Cy) in equation (16) can exclusively be 
expressed in terms of the nodal displacements and angles 
considering the requirements of continuity at the nodal points 


in x== 0) = Wy = Cy 
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OW, = = 
55 (= 9) Cee re 
1 
aa = — i GE We Urs Ee - 4 
w(x=a) = Wy = 1 9 ca cya 
Bld (5-4) = W = ¢c, + 2c,.a + 3c an 
aX Xo 2 3 4 


The above requirements can be expressed in the matrix form, 


q=-re (17) 
where 
A ee «() 0 
Oi 30 0 
= = 5 3 (17a) 
1 aoa a 
i) De 2a 3a) 
At ey 
Cc Wy: 
; Z a: 
Cs z lel) and a0 Sages Go] e) 
3 2 
lc, | x, | 


From equation (17), C can be expressed 


catlg (18) 


~~ 


me eehe new matrix G is defined such that 


Ger? (19) 


F tee od 


then, 


tQ 
II 


ome 20) 
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where 
1 0 0 0 
0 J. 0 0 
G= , 5 (25) 
a -3/a -2/a 3/a -lfa 
3 


Wit ayes” ya? 
The second derivative of deflection function becomes, 
a WK. 2c. + 6¢C,X 


2 3 4 


maach is in matrix form 


where H is the matrix defined by 





een, 
a*sw) (a? 7 7 oT 
| 2 | : 3 sao = OS) Ja (= Op ae Uo) Oy I = 00 LO, 
dx ax 


but, using the equation (20), 


2 2. 
|: a bed = (6q")(G> HY HG q) (24) 
dx ax sr 





The matrix expression for the first derivative of the 
deflection w can be expressed in the same manner, 

dw _ | 2 

ax = C5 ae 203X + 3C4x 


Hsing F and C 


RI 


= FC (25) 
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where 
2 


ao Cie 12x. 3x) (26) 
fealily, using equation (20) 
az = (6g )(G EF) ai) 


and 


G q) (28) 


From equations (24), (27) and (28), the Minimum energy 


equation (15) can be written in the matrix form, 


~~ 


a 
tbe f (8q") (G FE. F G q)dx 
0 ~_ —— wa 
a 
-be f (6q') (G? F’)dx = 0 (29) 
ee ie 


Noticing that the matrix gq, ge G, and G is not the 
Bpumection of x, but only of the function of dimension of the 
plate, and displacement and angle at the nodal points; there- 


fore, equation (29) can be simplified to 


a 
Sq: {EI ee uy (ut H)dx]G q 
0 oa 


T 


T a 
the G [/ (F F)dx]G q 
; q 


a 
= ae. ce [Sf (F") dx] } = 0 (30) 
0 





Pay 
in which the only integrations necessary to be carried out 


are the bracketed quantities, which become after integration, 


H)dx = (31) 


No gO Gac 12a%f 


Ae 0 1 2a 3a’ 
ee dx = 5 5 (eZ) 
0 O 2a 4a 6a 
Lo ane ea 9a | 
0 
he il 
S(F°)dx = (33) 
0 Desh 
[3a] 


Equation (30) is the final equation for one element 
formulated by the finite element displacement method. The 
First quantity is defined as a stiffness matrix for an 


aT the second is the additional stiffness matrix due 


element, 
to welding, and the third is the equivalent load due to weld- 


ing. They are as follows: 





Seatfness matxix for an element, 


a 
k= EL G ([/(H H)dxlg 
0 


2 Z 


ey es eh) ey 


6/a’ 4/a “6/Eo 27a 
= EI 


ya Siar Bye ~6/a* 


Ds 


p/n) ee Ge Ae 


Additional stiffness matrix due to welding, 


~ 
tl 


a 
be Gi) [ f(g! F)dx]G 
0 


0 0 0 QO 
Ore. 0 
oF Or =0: 0 


lo 0 0 be} 


Load matrix due to welding, 


a 
@e— be ca F ax] 
0 
0 
0 
0 
hee, 
| O 


28 


(34) 


(35) 


(36) 


At this stage, assembling of equation (30) has to be 


considered when the one-Span beam is divided by many finite 


elements, but the one-span beam element of Figure 4a is 


enosen to be two finite elements as in Figure 4b; then, only 


element 1 can be used to predict the deflection and angular 
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change at the nodal points by symmetric argument. In this 
Case, the simultaneous equations to be solved are (from 


equation (30)), 


6EL ,2 2 
éw, [—= (= w, + w -=—w.+w.)j] = 0 
1 ar a 1 x4 a 2 Xo 
ZEL 3 3 
éw —— (— w, + 2w -—w.+w)j = 0 
xy a Se xy a 2 Xo 
6EIL 2 2 
6w. [—~ (- —w.,-w + —w.-w.)j)] = 0 
2 5? a lk x4 a 2 Xo 
Ee. 3 3 abc = 
ow. a es Wal + we - Ps Wo + (2 + 5ET) Wx - bee.) ] — (6) 
2 ih 2 
(37a-d) 


From the simply supported boundary condition at node 2, 
Wo ="Q; therefore, OW. = 0. 


From symmetry condition at node l, Wi, a O; therefore, 
ale 


Furthermore, the displacement at node 1 and the angle at 
node 2 is not known, which means that the variations of these 
variables are not zero; therefore, the only way to satisfy 
equation (37a) and (37d) is the bracketed value in each equa- 


tion has to be zero, which leads to final equation for w 


and Wes ; 
2 
2 _ 
ai” Vice 
3: abc be _ 
ed (2 + ZBI "x, = 0) (Jean) 
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From equation (38a-b), displacement at the middle point and 


emgulac change at the edge expressed as, 


Ee -ado 
Le adi 1 ya) (39) 
Ce «CD + Bi/abc) (40) 


Y, 
It is worth to compare the results of equations (6) and 
(7) with equations (39) and (40). In equations (39) and (40), 


Se 
sw eee eee oe 


and a = 2/2 is to be used. Then, 


— 20 


mie? A4(L + 26I/ kc) (41) 


= Fo 
Bes | bt 2EI/te sete) 
where I is defined by unit width basis, ‘17° which is 


I= ¢°/l2. 
Equation (42) agrees with equation (7), as also does 

Begation (41), if x = 2/2 and @ of equation (7) is plugged 
into equation (1). The negative sign of equation (41) means 
that the coordinate and the positive angle at the edge has to 
be chosen, aS in Figure 4a. 

As shown, the results from the finite element method is 
exactly the same expression proposed by Masubuchi, Bese te": 
for the one dimensional deformation, and in this way, the 


equivalent loads and the stiffness matrix due to welding 


become easily visualized. 
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B. Two Dimensional Case 


Beeeeormulation of the Problem 

In the two dimensional analysis, all procedures are the 
Same as that of the one dimensional analysis except that the 
integrations are carried out by the numerical methods and the 
methods of assembling are mentioned. The deflection may be 
approximated by the bi-cubic function for one Bikawene , 
as in Figure 5. 

w= Cy + Cox ou C2yY + CAxY + c.x* + cey* + c4x" 


Z 2 3 3 Za 
+ Cox y + Coxy” + Cypy + CyyX ¥ + Cyox y 


3 Bane As a5 
+ Cy, 3xY + Cyyx y+ Cyex y + Ci gx Y (43a) 


fz 


2 Z 
Co = CyY sp 2C,X +5 3C5X “ 2C oxy os CoY 


ro) 
x 


2 a 3 Dimee. 
i 3C, 4% vet 204 XY alr C1 3Y ats 35 4X Vy 


2 Boas 
+ 2c, XY + 3c, (x y (43b) 


ma Z. Z 
C, + Cyx + 2EE Ya Cox” + ZO ex “at 3Cy oY 


6 9 


ajar 
cle 
[ 


3 
3 Z Z S 
+- C11 ate 2C, 5X Viet 3C4 XY ate 26 4% 4 


Saye 


252 
ap 3C) 5x We 3C, 6x y (43c) 


= 2 
Bey = C4 +3 2C 9X + 2C oY 4° 3¢54% + 4c, 5*Y 


Zane 


z. 2 2 
+ 3c,3y° + 6c,,x"y + 6c, ,_xy” + 9c, xy (43d) 








Lis = We 

fr = US % 4 

Lis = Wy 
y Lie = UT x44 


DIRECTION oF 
NODAL NUMBERING 


(0,0) 


q) = our 

ee ux 1 
Qa = Way 
14 = Wry 





on 


14 = Ws 
Lio = Wx 3 
eae Wy 3 


Vi22 Wxy 3 
(Q,b) 


Figure 5 Local Coordinate of a Two dimensional Rectangular 
Finite Element and Generalized Displacement Numbering 
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From equation (43a-d), the matrix equation can be defined 
such that at each nodal point, the generalized displacements 
(deflections, angles in x and y, and twisting angles) are 
expressed in terms of the geometry of the element and the some 
unknown coefficients Cy which is the form of, 

q=2TC (44) 
Weeee 1 1s 16 by 16 matrix, C is 16 by 1 matrix, and gis 
16 by 1 matrix defined as in Appendix A. 


meem equation (44), the unknown coefficient matrix, C€, 


can be expressed in terms of geometry and the generalized 


displacements, 

-- So (45) 
where 

G=r" (46) 


mewever, the actual calculation of the inverse of T is 
carried out by the computer in this analysis. 


(8,9,10) the strain energy 


From the elastic plate theory, 
stored in the plate during the elastic deformation is ex- 


meessed by 


a bil.2 Ne YD. 2 2 g 
U, = Ef fj + SB) = ati-y) |228- 28 - aay axdy 
Oe Ol ipox oy ox OY 


i) 


Meas can be written in the matrix form, Winicla ls 
7 72% of 
mee= =-—-f f K D K dydx (47) 


where 








34 





i= o- gw = =} c (48) 


9 2w 


easy 


in which H 1s given by Appendix A, 


‘eo 0 
FE | 
Dy = ee Vv 1 8) = B e D (49) 
12 (1-v~) 
1-v 
6 0 2 
in which 
7 Et? 
B= 2 
12(1-v~) 
and * 
ie 0 


© 
© 


LY 
2 


he, 


Therefore, equation (48) becomes after taking the 


Variation with respect to generalized displacement, 


ab i 
ee = f/f f (6K )8 DK dydx 
[oa 00' 
ab 
=f f (8C")(H’ 8 DH C)dydx 
00 
T m ab 
= 6g 8G ly su DH dyaxlG gq (50) 
0 0 
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The welding energy for the elements shown in Figure 6b 


can be written as, 


acy 2 a 2 
on =: é 5 (9,78, dx + : > (0 -@) dy 
Therefore, 
ey, el 2 c, b 2 
ey _ ow x _ ow 
Uy a, d We ay? eB ss BD: : Oe 3) a0 (51) 
which becomes, 
Te = Cy Se Cy. b OW A 
eee 2 2 sy) OP as Ay) 


1 2 
4 ee + C185] (1 
mad taking the variation, 


— ’ 
_ d0w| | OW dOW] | OW 
“9 c, : | | a anne: / | ee [Be | 


* baw salon 
2 Ss / age ax + C8. ‘ ae dy \ 53) 


in which the first bracketed term is the additional stiffness 
matrix due to welding, and the second is equivalent load 
matrix for the welding. In matrix form, equation (53) can be 


written, if some matrices defined as below: 


(54a) 
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——=F C=C F (54b) 


where the subscripts x and y denote that the integration 
should be carried out along the x- and y-direction, respec- 
tively, and the matrix F, and F., are defined in Appendix A. 


Therefore, 


d6w|{|ow ak Ak Jered Lars 
——||=—| = (dC Ee eC) = 6 G FO FG 55a 
| pe | [ae oy ey Ry ©) 212 aay 2 
dow] iow] _ ls él _ ee 
[2p (Se = (cC Je eee) = og (Gr ae (55b) 


Using equations (54a-b) and (55a-b), the welding term 


Cam be written as: 


a 
2 a T 
8, = oe CG - (FL FL)dxlG q 


b 
+ ogc, GIS (EY Ey)dyig g 


ne ~ 


a 
ac i T 
eeca C8 ah (Fo) dx] 


b 
= éq” ©. ol (EY) ay] (56) 


As in the one dimensional case already shown, the total 
Variation should be zero to be in equilibrium of this system. 


Therefore, 


6U,, = 6U_ + 6U.. = 0 
Ji p W 


wameh is from equations (50) and (56), 








a 


ab a 
Ab 
og") 4S" EB el Hq” D He dydiaa scam pooR dx 
oe 0 0 Y 0 ee 


b 
lt ab Ak 
pee s B fave g — G le toe si eee. 
0 
b om : 
+ C,8, : Bo ay | = 0 (Say) 


These are the 16 by 16 simultaneous eguations which have 
to be solved if only one element is concerned as shown in the 
one dimensional case, with the proper boundary and symmetry 
Semattions. From equation (57), the stiffness ke qeneaie< AF ONG &. 
plate element, | 


k= G 8 i i nu’ pH dydx G | (58a) 
0 0 


The additional stiffness matrix due to welding along the 


x-direction of the plate, 
T S At 
=e a she (c (58b) 
~WX ae ar 


The additional stiffness matrix due to welding along the 


y-direction of the plate, 


b 
=¢ C. FoF dy G (58c) 


k 
~Wy SY sane, 
The equivalent load matrix due to welding along the 


me-direction of the plate, 


a 
Q.= GC ens ye dx | (58d) 
~x oS MS) 0 


The equivalent load matrix due to welding along the 


y-direction of the plate, 


b 
oe CG C8, ; os dy 
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Figure 6 Various Finite Element Having Welded Edge 
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2. Assembling Procedure 
For all the finite elements except those having the 

welded edge, the stiffness matrix necessary to be assembled 
is the expression given by equation (58b), and for the finite 
elements having welded edges, three cases are caver weal. The 
first case 1S elements having welded edges along the x-direction 
as in Figure 6b; the seeone case is along the y-direction, as 
in Figure 6c; and the thira2 case is along x- and y-directions 
aS in Figure 6b. The assembiisuny procedure is the expression 
of equation (57) for the overall plate which is to be solved, 


which is described in FEABL user's AWA Tee 


Also, from 

the symmetry condition, which-will be described in the follow- 
ing section, only one-quarter part of the plate is actually 
being considered. 


Example of the assembling procedures for the two finite 


elements is given in Appendix B. 


3. Symmetry and Boundary Conditions 

After assembling, relevant symmetry and boundary 
conditions have to be applied to solve the simultaneous 
equations, as shown in the one dimensional case. Along the 
welded edge of x and y, the simply supported edge condition 
requires that all the displacements at that nodal point have 
to be zero, and for the Pereeran edge, the angular changes 
in the x-direction are not allowed, In just the same way, 
for y-directed edge, the angular changes in y are not allowed. 


The symmetry condition requires that the first derivative 
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4l 
along the symmetric axis should be zero. For example, the 
symmetry and boundary conditions for the case given by 
Figure 7 are as follows: 

Boundary conditions: 
Re node 6, 12, 18, 24, 30, 36, 42764 7eoe 
w= 0 
w. = 0 


x 
At node 40gayi. 57" 2, 53, 54 


Symmetry conditions: 


Memode 1. 2 ena 5 6 


w = QO 

; Y 
Memnodewlae  i-mero, WO, 25: S31. Spas 
sale 0 


Wie noae 1,16, 649 


Way = Q. 


4. Numerical Integration Method 

Usually in the finite element analysis, the integrations 
involved are carried out by the numerical methods, which can 
be easily calculated using the computer. For the present 
analysis, the nine Gaussian points are chosen for the rectan- 


gular finite avenent, - aao 


as shown in Figure 8c. In actual 
integration, it is very convenient to use the normalized coor- 
dinate system, as will be described below. From Figure 8a and 


8b, normalized coordinate system can be defined as 








42 


x = $(x + 1) (59a) 

y = Biy + 1) (59h) 
and 

dx = 5 dx | | (60a) 

ay = 2 ay (60b) 


Therefore, the arbitrary integration 


a 
ira ,b) = f P (x,y) dxdy (61) 
0 


(2 eae Ol 


becomes in the normalized coordinates, 


el. ae * 
fFla,b) = 22 ps gs P[S(X+) , >(¥+1)] (62) 


oe 1) 


meguation (62) can be integrated by numerically using 


9-GausSian points, which become 


b 3 3 
fla,b) = — EE WW. P(x,y) (63) 
i=l j=l J 
where 
=~ & 
x = 5 (Ks + 1) (63a) 
mee iy. + 1) (63b) 
2°35 
W.: W, = Ve Miley 8/9, Me 5/9 (63c) 
W, WS 2/97 87 a ee (63d) 
X.: X, = -Y15/5, X, = 0, X35 = "AIC AS) (63e) 
Y Y, = -/15/5, Y. = 0, Y, = Y1I5/5 (63f£) 
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Figure 8 Ordinally and Normalized Coordinates for an 
Element, and 9 Gausian Points for 
Numerical Integration 





44 


Therefore, using the concept of normalizing and 
numerical integration described above, the integrations 
involved in the present analysis such as equation (58a-e) 
can easily be written, which are programmed in Appendix D. 

From equation (58a) 


3 3 
2 Gn: WW. P(x,y)]G (64) 


k= ( 
i=l j=l J 


where 


p(x,y) =H’ DH 


~ 


x, y, Ws, Was X. and oe are defined by equation (63a-f). 


neem equation (58b), 


aC 7 3 
Pea (OS LOW, PL Ge ydiE (65) 
i=l] 
where 
Ae 
P.. (x,y) = ae Fy 
ass 
x = aX + 1) 
‘ee 


Ws and Xs are defined by equation (63c) and (63e), respectively. 


From equation (58c), 


bG 3 


= = ¢*I Wigs 2) leap) ue (66) 
a al evry ~ 
J 
where 
T 

pe (x =e he abe 

zy! ry) Be, by 

m= a 


kK 
Tl 
NO 
K 
J. 
+ 
i 








45 


W. and Ke are defined by equations (63d) and (63f), 
respectively. 


From equation (58d), 


adic T 3 | 
ee Gl Se Gy) JG (67) 
i=1 
where 
T 
Pox %¥) eX 
et 
—— 5 (Xs + 1) 
ca b 


W, and X, are defined by equation (63c) and (63 ), 


mespectively. 


From equation (58e), 





boc, rr 3 
— ete Ws Po (x,¥)IG (68) 
where 
Alt 
P ae F 
ae y) a 
xe a 
b 
=>(¥. +1 
y 5 | - ) 


W. and “ are defined by equations (63d) and (63f). 
As noticed above, numerical integration using the 
normalized coordinate can simplified the computing procedure 


Byetactoring out a and b. 
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me Computer Programs 

Computations of the ordinally, additional welding, and 
equivalent welding load matrix of equation (58a-e) have been 
programmed by using the numerical integration technique with 
normalized coordinates aS listed in Appendix D. In addition 
to these, FEABL Braoreresmas in assembling of the relevant 
matrices and solving the simultaneous equations with proper 
boundary and symmetry conditions are also used and listed in 
Appendix D. 

The Programs listed in Appendix D are only relevant to 
the panel structure welding with rectangular finite elements. 
Further, the nodal numbering system should follow the way 
given by the sample analysis, which means that the numbering 
should start from origin to y-direction, as shown in Figure 7. 

Input data for this analysis are the overall dimensions 
of one panel structure in the x- and y-direction, thickness 
of the plate, Young's modulus of elasticity, Poisson's ratio, 
number of elements divided in x- and y-direction, equivalent 
welding moment in x- and y-direction, and free joint angular 
Change at the edge. The input data cards are described in 
Appendix C. 

Outputs are the displacements, angular changes in x- and 
y-directions, and the twisting angles at each nodal point, as 
shown in Appendix E. 

One thing to be noticed here is that for a quarter part 
of the plate, as shown in Figure 9, the proper boundary and 


Symmetry conditions at the nodal points are to be generated 
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by subroutine HOLD given in Appendix D; therefore, as long as 
boundary conditions of simply supported and symmetry conditions 
as described are concerned, it is not necessary to be consi- 
dered in the input. Therefore, no matter how the number of 
finite elements are to be chosen, the only input data will be 
as described in equation (69a-h). Furthermore, to use this 
program, the coordinate and a quarter part of the plate should 
be chosen as the same way as shown in Figure 9. 

The units used in the sample program are pounds, inches, 
and radians. As results, for the output data, the deflections 
are in inches, the first derivatives (angular changes) are in 
radians, and the second derivatives with respect to x and y 
(twisting angles) are also in radians. 

-Program language used is FORTRAN IV, and the computations 


are carried out by an IBM 370/155 system. 
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IIL RESULTS 


As computer inputs, one steel panel structural dimension, 
as shown in Figure 10, for which the experiment conducted and 
weight of consumed welding rods are given, 3 gr/em., ‘7 '+?? 
are used to predict the deformations, angular changes, and 
twisting angles at the nodal points. Here, the value of the 
One dimensional free joint angular change along the welded 
edge, Ses has been read fz27:the results of Figure 3 and the 
equivalent constrained welding moment, C, is being evaluated 
uSing equation (8). The numerical values of the input data 
mee as tollows: 

Half-length of the plate in the x- and y-directions, 
respectively, L, = 16 inches and L, = 10 inches. 

Number of finite elements in the x- and y-directions, 


respectively, M = 8 and N = 5. 


Equivalent constrained welding moment from equation (8) 


where W is given by 3 gr/cm from experiment. 


Therefore, 


<= SONS Wega eae 


Petsson's ratio: 
ba- 0.3 
Young's modulus: 


6 


Ee 30 x 10° psi. 
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Thickness of the plate: 


i= 95/8 anehy. (69 a-h) 


The results of the analysis using the above data are 
listed in Appendix E, and these results are plotted in the 
Figures 11 and 12 to be compared with experimental results. 
me Can be noticed from Figures 11 and 12, the finite element 
results are lower value than the experiments, and the maximum 
deflection at the midpoint of the plate which is the most 
important deflection in reality is only a half of that from 
the experiment. 

To visualize the behavior of the deformations as 
mnccions of ae and C, Series of computations are conducted 
with different combinations of oe and € as shown in Tablewe 
Yet remaining, the rest of the input data are unchanged, and 
the maximum deflections at nodal point 1 are plotted in terms 
of ©, and C as in Figures 13 and 14, and the maximum angular 
Changes at node 6 and 49 are plotted in the Figures 15 and 16. 
Only the maximum deflections at node l, and the maximum 
angular changes at node 6 and 49 from the results of series 
computations are tabulated in Tables 2 and 3. 


3 and 


Finally, computed results using ee = (55) < ale 
C= 4.2 x Or are listed in Appendix E, and also plotted in 
Figures 11 and 12 to compare with the experiment. As can be 
seen from Figures ll and 12, these results are very close to 


the experiment and the reason of choosing the above 8, and Cc 


werd be discussed in the following section. 
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iV DISCUSSIONS OLS a sUnr. 


As has been mentioned in the previous section, the 
results using C = 1.13 x 10° and @ = 45 x 10°° are only a 
half of the experimentally measured values as shown in 
Figures ll and 12. The reasons of this difference are not 
yet known. These may result from the assumption of the elas- 
tic behavior of the welded deformation, or from the uncer- 
tainty of the experimental measurement itself because the 
deflections and angular changes which are only of the order 


inches and La radians, respectively. 


of magnitude of 10° 
met, Other possibilities are from the value of 86 havc, 1G, 

To investigate the behavior of the deformation with 
respect to 86 and C, series of computations have been con- 
ducted. The plotting of maximum deflections at node l 
versus C with changing 0, as in Figure 13 shows some insight 
Of deflection behavior with varying C. In here, it is ob- 
vious that the deflections are more rapidly changed with 
changing =, than C changes, which means that the deflections 
are more strongly dependent function of oe rather than C. 
Furthermore, the deflection reaches rapidly to an asymptotic 
value as C increases, which can be seen clearly from the one 
dimensional deflection, equation (41) where the term 2EI/ &c 
becomes small compared to 1 as C increases; therefore, in the 
- limit case, the deflection can be expressed only by 


w= £0 /4. 








62 


This behavior of the deflection with changing C are 
clearly shown in Figure 13. Also from Figure 14, in which the 
maximum deflections are plotted with respect to ee it is 
obvious that the deflections are rapidly increased as as 
increases. | 

If assumed that the difference between experiment and 
computed results are caused from the possibility of different 
values of i and C, then the question is how to choose Oe and 
C, such that the computed maximum deflection becomes the same 
to that from the experiment. From Figures 12 and 13, the com- 
binations of ne and C to satisfy the given value of maximum 
deflection from the experiment can be predicted. However, 
the best combination is not known, but for the known experi- 
ment value of maximum deflection it may be safe to choose oh 
and C such that the changes of both - and C from the cal- 
culated value be minimized. In other words, superposing 
Figures 13 and 14 as shown in Figure 17, the point which 
gives the same reading of A and C in both ways as shown by 
the small circle, which is read as oe = 55 x O02 and 
C= 4.2 x 18, will be the optimal combination to satisfy 
the given maximum deflection. The comparison of the results 
using this 8 and C to the experiments, as in Figures 11 and 
12, shows that this approach gives reasonably good results. 

As has been seen from above, there exists many 
uncertainties for the two dimensional welding analysis, but 
the welding deformation phenomena can be predicted with a 


reasonable degree of accuracy by the finite element analysis. 


eae 
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Furthermore, the empirical formula of C given by equation (8) 
may not be simply applied to the two dimensional case, and 
also does the value of O5° If it is so, meaSuring only the 
maximum deflection in future experiments, the best combination 


of a and C can be predicted as discussed. 
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V CONCLUSIONS AND RECOMMENDATIONS 


mo Conclusions 


From the discussion of the previous section, it can be 


concluded as: 


de. 


The finite element analysis for the two dimensional 
welding deformations which are assumed to be elas- 
tic deformation gives reasonable prediction of 
deflections. 

The searching technique of optimal combination of 

OG and C for a given deflection can be used as shown 
already, if there exists the uncertainty of 0, and C. 
fimeche future experiments, this analycois Giveseene 
order of magnitude feeling of the deflection mea- 
sured after welding. It is extremely difficult to 
measure correct deflections‘in experiments because 
the order of magnitude to be measured is very small; 
therefore, the inherent uncertainties always exist. 
The magnitude of importance of the welding deformation 
variables such as a and C are clearly seen, there- 
fore, more care should be concentrated in measuring 
oP rather than C. 

Simple empirical formula for C given by equation (8) 
may not cover the wide range of the values of C; 
furthermore, as can be derived from equation (41) it 
may be logical to express C in terms of maximum de- 


flection, free joint angular change, and rigidity, 
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such as: 


C= eee - 2EI/2 (70) 


in which De is a function of a given welding 
Send1tleon and Wess is the constrained maximum deflec- 
GLOn at midpoint o£ the span for the Gamleseoncicion 
of the welding as of OO" Further, the elastic 


rigidity and plate dimension terms are expressed by 


Ene term Of Zee. 


B. Recommendations 


In this analysis, the assumption used is the elastic 
Mehpavior of deformations. As stated in the conclusion, it is 
still believed reasonable to use the elastic theory, but for 
more accurate prediction, local plastic deformation phenomena 
should be included in the analysis. 

If the elastic assumption is to be used, more careful 
analysis to evaluate the value of oR and C are recommended in 
future research. Also, a more rational approach to express 
C is recommended. 

In future research, it is recommended that series of 
experiments be conducted to measure the maximum plate defor- 
Mation only, and express C and 8, in terms of welding and 
plate dimension. In this way, the computer inputs are always 
available for any prediction of deformation for the two 


dimensional panel structures. 
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APPEND I Xe 


Example of Assembling 


As an example, the assembling procedure for a quarter 
part of the plate divided by two finite element Snow ae 
Figure 18a will be discussed. 

For the element 1, as shown in Figure 18b, the stiffness 
matrix in a local coordinate can be written from equation 
f5ca, Db), 


ioe > kK + ik (A-1) 
me 


where KL wsol6 by 16 symmetric matrix. And fem 
1 
equation (58d), the load matrix for the element 1 is, 


Q., 7 & (A~2) 

where Q, is 16 by 1 matrix. By the same way for element 2, 
Ib 

the stiffness and load matrix can be written using equation 


(58a, b, c) and (58d, e), respectively, 


k, =k+k vk (A=3) 
ie = oe a By (A-4) 


where Kk, is 16 by 16 symmetric matrix and Q. vs, lo by st 
5 = 


2 
matrix. 


At every nodal point there are four unknowns (w, Bh Tat 


Wary)? therefore, after assembling element 1 and 2, the 
equations have to be solved with proper boundary and symmetry 
conditions are 24 by 24 simultaneous equations, which is in 


the matrix form: 





a. 


6g [Ka - Ql = 0 (A-5) 


where the matrix K is 24 by 24 symmetric matrix, Q is 245Dyo4- 
q is 24 by 1, and i. is 1 by 24 matrix. 

All the quantities in equation (A-5) are defined by the 
global coOrdinate, therefore, the assembling procedure is to 
i@edee the quantities defined by equation (A-l, A=-2, A-3, A-4) 
to the proper position in equation (A-5). This operation can 


easily be established comparing each nodal points numbered in 


local coordinates to that of global coordinates: 


Element 1, nodal points: Element 2, nodal points 
Coordinates -Coordinates 
Hocal global local global 
1 1 1 3 
2 5 2 5 
3) 4 5 6 
4 2 4 4 


And noticing that numbering of the generalized displacement, 


Qs follows the relation, 

i=4m-g (A-6) 

j= 4n-h (A-7) 
where i and j represent the index for generalized displacements, 
mand n for nodal numbers, g and h are set to be 3 for w, 


Zz for wl, 1. £6nsw angen con. yee 
x yY xy 


For example, one of the stiffness matrix elements at 


nodal point 4 in global coordinate can be calculated. 
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L ae 
"SE on a 


a. Nodal points numbering in global coordinate 





XK 


b. Nodal points numbering in local coordinate 


Figure 18 Elements and Nodal Points Numbering in 
Global and Local Coordinate for Sample Assembling 
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PaOmrcouatlione (A ©) ano mh) 
fo 16 - ¢ (A-8) 
7 = 16 -h (c=) 
But the same node 4 becomes node 3 at element 1 and node 4 
at element 2 in local coordinate. 
For element 1, 


(1) 4 eZ ee (A-10) 


ll 


oP = 12 - ib (A-11) 


and for element 2, 


(i), = 16-g (A~12) 


i 


(3) 5 16 -h (A-13) 


Therefore, if g = 1 and h = 2, then from equation (A-8), 


(A-9), and (A-5), 


ey «CS 15,14 (A~14) 
From equations (A-10), (A-11), and (A-1) 
Sa (A-15) 
i, 10 
ema from equations (A-12), (A-13), and (A-3), 
ss (A-16) 
15,14 
Therefore, element of the stiffness matrix in global 
coordinate becomes, . 
Kis qa = Ke, eet Ke, (A-17) 


IE 5 Law, 57 





All other stiffness matrix element can be expressed by the 


same manner described above, and the programs of assembling 


the stiffness matrix are available in FEaBL. ‘2? 
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APPENDIX 


DeSCEBpELON of Input Data 


Input data for the programs listed in Appendix D are: 


% ~— Oo. , «OUD, h6Ut, E, Loy i M and N as described in 


ae O ve 
section III. In addition to these, control card for calcula- 
tion of different cases in each job with changing above input 
variables are considered. By using this calculation control 
input card in each job, series of calculation with changing 
Some of the input variables can be calculated, and the defor- 
mation behavior with respect to these changing variables can 
easily be visualized by plotting the deformation versus these 
variables. 

Detailed input data card arrangement and format are 


. 


imested in Table 4. 





Table 4 


inpuc Data Layout-formm 


Card Pipues rogram 
Column Format Symbol 
Card Number 1 

1-5 iS NCASE 
Card Number 2 
moe B15. 7 LX 

Moeeo) b15.7 LY 

31-35 15 M 

36-40 I15 N 
Card Number 3 

enor 2.7 Cx 
mo—-24 B1l2.7 CY 

25-36 E12.7 TESS UID AN ve 

37-48 El12.7 GNU 

49-60 E12.7 THK 

61-72 E12.7 E 


Definition. and Unit 


Number of cases involved 
in each job 


Length of 1/4 plate in 
x, inches 


Length of 1/4 plate in y 
Number of elements in x 


Number of elements in y 


Equivalent constrained 
welding moment as in 
Figure 6, lb-in/in rad. 


Equivalent constrained 
welding moment as in 
Figure 6, lb-in/in rad. 


Free joint angular 
Change, 8,5, radians 


Poisson's ratio, D 


Plate thickness, t, 
inches 


Young's modulus, E, 
lb/in.2 
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Sample 
_Input 


14 


16. 0n4200 


10.0E+00 
8 
5 


4.20E+04 


4.20E+04 


SO OB Ue 


OF SUE00 
Ono7 Sa UW 


30. CETOG 
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